WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau llON 



PCT 

INTERNATIONAL APPUCATIQN, PUBLISHED UNDER THE PATENT COOPERATION TREA TY (PCn 

WO 96/25673 

22 August 1996 (22.08.96) 


(51) International Patent Classification 6 : 
Cr S ' 5/14, 1/02, 11/02 


Al 


(11) International Publication Number: 
(43) International Publication Date: 


(21) international Application Number: PCT/GB96Y00270 

(22) International Filing Date: 6 February 1996 (06.02.96) 


(30) Priority Data: 

08/389,263 


16 February 1995 (16.02.95) US 


(71) Applicant: NEXUS 1994 LIMITED [GB/GB]; 7-10 Chandos 
Street, London W1M 9DE (GB). 

(72) Inventors: YOKEV, Hanoch; 8 Hazanchanim Avenue, 52341 
Ramat-Gan (IL). PELEG, Shimon; 11 Anchilevich Street, 
45285 Hod Hasharon (IL). MEIMAN. Yehouda; 31 Ben 
Eliezct Street, 75299 Rishdn Letzian (IL). PORAT, Boar, 
39a Shimshon Street, 34678 Haifa (IL). 

(74) Agents: BERESFORD, Keith, Denis, Lewis et al.; Beresford 
& Co.. 2-5 Warwick Court, High Holbom, London WC1R 
5DJ (GB). 


(81) I*signated States: JP, KR, SG. European patent (AT, BE CH. 

• D *> ES, FR, GB, GR. IE, ITLU, MC. NL, h^E) 

Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. J { 


CO 


o 

o 


100 


(54) Title: DIFFERENTIAL RANGING FOR A FREQUENCY-HOPPED REMOTE POSITION DETERMINATION SYSTEM 
(57) Abstract 

A differential ranging location system is de-fj 
scribed which uses a modified time-of-arrival tech- ijf 
nique to determine the location of a frequency- ,t! 
hopped spread spectrum radio signal. The trans- i 
mitter simultaneously transmits two radio frequency i 
carriers h aying d ifferent nxquerSes"SuclT that a S 
phasedft ference is^b"seVve*r oerweeh'the two carri- 
ers at a distance from the transmitter. The phase dif- 
ference is proportional to the range from the trans- 
mitter that the carrier signals have observed. The 
two carrier signals from the single transmitter are 
received by at least three and in special cases four 
base stations which calculate the differential time 
of arrival based on the phase differences of the re- 
ceived carriers. The calculated phase differences 
are then sent to a central location which locates the 
position of the transmitter based upon a planer hy- 
perbolic location algorithm. Since a large number 
of narrow band frequencies across a wide spectrum 
are used in a frequency hopping transmission pro- 
tocol, the accuracy in determining the phase differ- 
ences between the two carriers is increased and the 
immunity to interference is also increased. 
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DIFFERENTIAL RANGING FOR A FREQUENCY-HOPPED 
REMOTE POSITION DETERMINATION SYSTEM 

FIELD OF THE INVENTION 

The present patent application relates to radio communication 
and transmitter location system. 


10 


BACKGROUND OF THF 1N1vpm T1 ^ n 
Location of radio transmitters is a technology that has been 
used for many years and in many applications. For radio location, the use of 
mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio signal 
by triangulation of the radio signals. Triangulation is a technique where a 
plurality of fixed base station receivers each determine the approximate 
15 direction of the radio signal from the transmitter. The approximate direction 
lines are then drawn on a map and the location of the transmitter is defined at 
the cross point of these lines. Since the direction lines are approximate, and 
typically three fixed base station receivers are used, the intersecting area 
between the three lines is usually shaped like a triangle since the lines rarely 
20 meet at a single point. The transmitter is then located within the triangle. 

Triangulation of radio signals is a costly undertaking since at 
least two receiving antenna stations must be situated at different diverse 
geographic locations to perform even rudimentary triangulation. Effective 
triangulation is more realistically achieved with three receiving locations and 

25 accuracy is increased even more with four receiving towers. The more 
receiving stations constructed, the better the accuracy in locating the 
transmitter, but also the great the cost. 

One such direction finding technique is described in US Patent 
No. 5.379,047 entitled "REMOTE POSITION DETERMINATION SYSTEM" 

30 issued January 3, 1 995. In that patent, the direction of the incoming radio 
signal is determined by using an array of receiving antennas to measure the 
phase difference between the antennas of the array and thereby calculate the 


WO 96/25673 


PCT/GB96/00270 


2 


incident direction of the radio carrier signal. Although this patent has many 
novel and useful advantages, it requires that each base station use an array of 
antennas and at least two and typically three bases. Quite often the cost of 
building and maintaining these base stations can be quite high, especially if 

5 they are located in urban areas where roof-top space must be leased and the 
cost of the lease is determined by the number of antennas placed thereon. 

Another disadvantage of triangulation for location of radio 
transmitters is the construction of such a system in sparsely-populated areas. 
In the example of a highway across a desert area, it is cost prohibitive to 

10 install a large number of base stations to track a narrow stretch of road. The 
base stations would have to be placed on both sides of the road and spaced 
quite far back from the road. 

There is a need in the art to locate a radio transmitter using at 
least three and in special cases four one-antenna receivers. A receiver with 
15 one antenna will minimize the construction and lease costs of a location 
system. There is a further need in the art to determine the location of low 
power, frequency hopped spread spectrum radio transmitter. 

SUMMARY OF T H E INVENTION 
20 The present invention solves the above-mentioned problems 

with locating transmitters and other problems which will be understood by 
those skilled in the art upon reading and understanding the present 
specification. The present invention uses a modified time -of-arrival technique 
to determine the location of a frequency hopped spread spectrum radio signal, 

25 The transmitter simultaneously transmits two radio frequency carriers having 
different frequencies such that a phase_djffere n ce is observed between the two 
carriersju a distance from the transmitter. The phase difference is *~ 
proportional to the range from the transmitter that the carrier signals are 
observed. The two carrier signals from the single transmitter are received by 

30 at least three and in special cases four base stations which calculate the 
differential time of arrival based on the phase differences of the received 
carriers. The calculated phase differences are then sent to a central location 
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which locates the position of the transmitter based upon a planer hyperbolic 
location algorithm. Since a large number of narrow band frequencies across a 
wide spectrum are used in a frequency hopping transmission protocol, the 
accuracy in determining the phase differences between the two carriers is 
increased and the immunity to interference is also increased. 

BRIEF DESC RIPTION OF THF QRAWINfiS; 

In the drawings, where like numerals refer to like components 

throughout the several views. 

Figure 1 is a diagram of an integrated frequency hopping spread 
spectrum communication and location system; 

Figure 2 is a block diagram of a reverse paging and mobile 
location system using the infrastructure of a standard paging system; 

Figure 3 describes the synchronization and message format of 
1 5 the outgoing paging signals from the base stations; 

Figure 4 describes the format of the frequency- hopped spread 
spectrum signal transmitted by the remote mobile units; 

Figure 5 is a hyperbolic plot of the differential ranging of the 
present location system; 

20 Figure 6 isa di *gram of the differential time-of-arrival aspects 

of the present location system; 

Figure 7 shows the computations used to perform location using 
differential time-of-arrival ranging; 

Figure 8 describes in more detail the principle of transmitting 
25 two simultaneous frequencies to two receiving base stations; 

Figure 9 is block diagram of the differential ranging transmitter: 

Figure 10 is a diagram of the IF circuit of the differential 
ranging receiver in the base station which uses two separate convenors: and 

Figure 1 1 is a phase plot showing the effect of sources of errors 
30 on the differential ranging technique. 
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DKTAIU-I) DHSCmPTION ()|: Till: l^ -n ^p KMMOniN/1WKn - 
In the following detailed description of the preferred 
embodiment, reference is made to the accompanying drawings which form a 
part hereof, and in which is shown hy way of illustration specific preferred 
embodiments in which the inventions may he practiced. These embodiments 
are described in sufficient detail to enable those skilled in the art to practice 
the invention, and it is to be understood that other embodiments may be 
utilized and that structural, logical and electrical changes may be made 
without departing from the SCODe of the presem invemions ^ 

following detailed description is, therefore, not to be taken in a limiting sense, 
and the scope of the present inventions is defined by the appended 
claims. 

System Overview 

In the preferred embodiment of the present invention, the base 

15 stations and central site _„ ■ .. „ ., , 

are similar to the base stations 

and central site described in U.S. Patent Number 5,430.759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 

PAGING SYSTEM" issued July 4, 1995 an H ih* r™~* • . 

j > yvD. ana the processing structure of this 

system operates similarly to the present system/ The present invention allows 
the accurate location of the frequency hopped; spread spectrum transmitters 
even in noisy and multi-path environments such as those found in urban 
environments and irregular terrain environments. 

The present invention is directed to a location system for use in 
locating remote mobile units (RJvlU's). The remote mobile units may be 
25 located in motor vehicles, located on the person of people, carried within 
containers or packages, or any number or variety of mobile carriers. The 
location system comprises a plurality of base stations (at least three and in 
special cases four) which can transmit signals to the remote mobile units and 
receive signals from the remote mobile units. The base stations may send 
messages and information to the remote mobile units, the base stations may 
send alert or interrogation commands to activate the remote mobile units, or 
the remote mobile units may be activated locally without any prompting. 


20 
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Thus the remoie mobile units may he activated locally (lor example by the 
wearer or vehicle) or remotely by the base station. The remoie mobile units 
may be used by a person as a reverse payer, as an emergency locator or as u 
communication device. The remoie mobile units may also be used as part of 
5 a vehicle to locate the vehicle if stolen, for vehicle tracking or as an integral 
communication device. The remote mobile unii may also be activated in the 
case of an accident by a deceleration switch. ABS system, air-bag or other ' 
crash-activated switch to automatically summon help 

The remote mobile units receive messages as a standard paging 
10 device over licensed airwaves using a standard paging infrastructure. The 
remote mobile units, when activated, transmit low-p OWer (less than one watt), 
frequency-hopped, spread-spectrum communication signals. The transmitted 
signals from the remote mobile units are received by the base stations which 
are then used to locate the remote mobile units. The base stations are 
15 equipped with a single receiving antenna. The location of the remote mobile 
units is determined by differential ranging using phase difference comparisons 
between two carrier signals simultaneously transmitted by the remote mobile 
unit. Unique algorithms are performed on the received signals at the base 
stations to eliminate the ambiguity caused by signal interference, phase 
20 ambiguity caused by errors and noise and multi-path reflections. 

The frequency hopping transmitter of the remote mobile unit 
produces two sine waves for each frequency hop. The phase difference 
between the two sine waves will vary with distance from the transmitter which 
is proportional to the range from the transmitter. There is an ambiguity factor 
25 in the range determination due to the phase differences between the two sine 
waves. The ambiguity is minimized by transmitting at least 53 hops per 
transmitted message so that there is a different phase difference within each 
frequency hop. Each hop change produces a step in frequency such that 53 
steps of various step distances in frequency are produced. 
30 At the "ceiver. the frequency hopping carriers are 

downconverted. sampled and digitized for processing in the digital domain. A 
range calculation is determined by comparing the phase differences of the two 
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sine waves of the two frequencies within each hop. The ambiguities are 
minimized by comparing the phase differences across the 53 hops in 
synchrony with the hopping protocol. The differential range between two 
bases is determined by collecting the phase differences from all the hops in 
two bases and performing in the central station an Inverse Fourier Transform 
(IFT) which will produce a single value lor the differential range between the 
two bases to the transmitter. 

Each base station is equipped with a very accurate rubidium 
atomic clock so that the exact time of arrival of the frequency hoppine signals 
is determined with reference to a very accurate standard time. Those skilled 
in the art will readily recognize that any highly-accurate clock means mav be 
substituted for the preferred rubidium clock such as a GPS receiver. The 
central station compares the time of arrivals with t..e knowledge that each 
base station has made its calculations based on the same time standard The 
time of arrival information is then used to plot possible locations of the 
transmitter using a hyperbolic plot. The comparison of the time of arrivals 
from two base stations will produce a hyperbolic plot where the base stations 
are the foci. By plotting the time of arrival numbers on the hyperbolic plot 
and crossing the hyperbolic plots between two or more base stations, the exact 
location of the transmitter can be determined 

To assist in calibratine the different;^ • 

"ufc. me airterential ranging system with the 

local base stations, a transmitter having a known location is used to 
periodically transmit a signal which is then used to synchronize the base 
stations. The calibration signals is used to adjust the local atomic clocks of 
the base stations. It is possible to use less accurate clocks in the present 
system by increasing the number of calibrations performed in a dav . In a 
calibrated system the accuracy of location has been measured at 26 meters at 
minimum received signal strength of -138dbm when operating within the 
parameters described below. 

The problems associated with multipath reflections and fading 
are known in the art. The multi-path problems and fading problems are 
usually frequency dependent so the use of a large number of IVequencv hops 
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vinually eliminates these problems. The annuity problems are also solved 
in a fashion s.milar to that described in the aforementioned eopendin,. paten, 
application. 

Summary of Hyperbolic Location wit hj^er cn tial Time of a.;,,,. 

The operation of the present differen^m^ system mav 
also be termed differential time-of-arrival (TOA) since the differences in the 
arrival limes at the base stations are used to determine the position of the 
receiver. As is described in more detail below, the TOA of the two carrier 
frequency signals for each remote mobile unit are estimated at each base 
station for each hop. The TOA is proportional to the range R of the signals. 

The operation of the present system is basically performed in 
two major steps. In the first step, estimated range K is determined for each 
base station BS n , where n is the base station number. Due to the inaccuracy 
15 of the clock in the transmitter, the range by itself Is not an accurate indicator 
of the distance of the remote mobile unit from the transmitter. Since the base 
stations are equipped with very accurate clocks which are periodically 
synchronized to a common time base, the base stations can report the exact 
time of arrival of the signals. 

20 In the second major step, the time of arrival information from 

each base station and the range information (calculated across the multiple 
hops using an Inverse Fourier Transform) are sent to a central station where a 
range differential between two base stations is plotted, as shown in Figure 5. 
The range differential forms a hyperbolic plot, f or example 501a and 501b 

25 between base station BS, and base station BS 2 . We know that the range R, is 
greater than the range R 2 so that the hyperbolic plot 501b in dashed lines in 
Figure 5 can be ignored. A second hyperbolic plot 502a and 502b is formed 
between base station BS 2 and base station BS,. The intersection point 
between the two hyperbolic plots is the location of the remote mobile unit 

30 RMU. 
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Comparison of IWernuin ) Kaimin» m CPs: 
The operation of the present dilTcreniial rani-ine system is quite 
ilit lcieni horn the GPS (Global Positioning Satellite) system. I he CiPS is a 
navigation system that provides loeation for the vehicle that carries it. In 
5 order to monitor the vehicles, the system requires a communication system for 
delivering the location to a central office. The navigation process of the GPS 
is virtually the inverse process of the dual frequency differential ranging 
system of the present invention. The theoretical accuracy of the present 
invention is only dependent on the bandwidth, signal to noise ratio and some 
10 ambiguities described below. Although commercial GPS systems are made 
inaccurate on purpose, a comparison can still be made. 

There are several facts in favor of the present differential 

ranging system for the task of remote mobile unit location over GPS. The 
GPS system requires direct line of sight to the satellites, since it has a very 
1 5 sensitive receiver. Differential ranging is transmitting in the 900 MHz band 
and thus does not require line of sight: it could operate from within buildings. 
Differential ranging does not require any other device other than the low cost 
RMU on the vehicle which makes it significantly lower cost. The power 
consumption of the differential ranging system is significantly lower since it 
20 requires only a low power transceiver, thus it is. operates on batteries for a 
long time (over 1 month on commercial batteries). The differential ranging 
system operates using frequency hopping and is less susceptible to multipath 
and interference than the GPS which operates in direct sequence spread 
spectrum. From measurements performed on the present direction finding 
25 system, it is by far less susceptible to ground multipath than the GPS. The 
complex receivers of the system are only in the base stations. The 

complex receivers of the GPS are in every vehicle 

Frequency Ho pping System nv».rvi»., 
30 In the preferred embodiment of the present invention. 

transmitters and base stations similar to those used in the present invention are 
described in the aforementioned copending U.S. Patent Number 5,430,759 
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entitled "I.OW-POWI-U FRI-Ol MiNCY-l lOPIM-n SPRKAD SI>1-(TI«IM 
RKVKRSK PAGINC, SYSTKM" issued July 4. 1005. The infrastructure of 
this reverse paging system serves as the basis ol' the present :U ui the 

reverse pagers of this system operate similarly in the present system using a 
5 single carrier frequency instead of the two carrier signals of the present 

invention. The system allows the accurate location of the reverse 

pagers even in noisy and. multi-path environments such as those found in 
urban environments and irregular terrain environments 

The environment in which the present system j s operated is 
10 shown in Figure 1. The system »s designed to operate in a 

communications environment known as NexNet™ in which a plurality of base 
stations service a particular area using the NexNet™ communications system. 
Each remote mobile unit receives standard paging messages as a downlink and 
can provide talk back messages as an uplink channel using frequency hopping 
1 5 spread spectrum radio frequency communication. The paging transmitter of 
the base stations 113 typically transmits with 200 Watts of power to send 
messages using a standard paging format such as POGSAG or other industry 
standard formats. Each remote mobile unit contains a paging receiver with a 
sensitivity of -1 15 dBm. a. frequency-hopping dual frequency transmitter with 
20 a transmit power of one Watt and double side band modulator: Each base 
station includes a receiving antenna for the frequency hopped carrier signals 
and special processors for determining the phase difference between the two 
carriers of the double side band modulation. 

Figure 2 depicts the major components of the two-way paging 
25 system in the aforementioned U.S. Patent Number 5,430,759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 
PAGING SYSTEM" issued July 4, 1995. The reverse paging terminal 110 at 
the central site operates to provide synchronization and messaging information 
through the paging terminal 1 16 to the reverse pagers 100 (also known as 
remote mobile units 100) via direct links to the base stations BS„ BS 2 and 
BS n through ground based radio links (not shown) or through a satellite 
uplink/downlink using a geostationary satellite (not shown). The base stations 
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BS,. BS ? mid BS, include iransmil and receive towers 113a. 113b and 1 13c. 
respectively and base station terminals 201a. 201b and 201c. respectively. 
Terminals 201a-201e are required for producing the accurate synchronization . 
information needed to be transmitted to the remote mobile units and for local 
5 processing of the received messages for transmitter location. This 

synchronization information is used to coordinate the frequency hopping 
transmissions and to coordinate the response of messaging from the plurality 
remote mobile units 100 so as to minimize collisions within groups of remote 
mobile units and eliminate collisions between groups of remote mobile units. 

10 

Base to Remote Synrhmn;-,,,;^ 
Standard paging messages sent from the base stations BS,. BS, 
and BS. to the plurality of remote mobile units 100 are. in the preferred 
embodiment, sent as digital data encoded in the POCSAG paging standard. 
15 These messages may be used to interrogate the remote mobile units to activate 
the. remote mobile unit to allow the base stations to begin the location process. 
. Typically the paging channel has a center frequency of 143.160 MHz, with an 

NRZ FSK data rate of 512 bps or. 1200 bps. Other bit rates such as 2400 
baud (bps) are also feasible. Figure 3 describes the POCSAG paging 
20 communications protocol as modified for use by the preferred embodiments of 
the present invention. In the top line of Figure 3. a greatly compressed time 
line of digital data transmitted according to the POCSAG protocol is shown. 
Batches or messages are transmitted in groups as shown in the details in the 
subsequent lines below the top line of Figure 3. In the second line of Figure 
25 3. a 1.0625 second interval (for 512 baud) is shown in which 544 bits are 
transmitted as a single batch. The batch is preceded by a synchronization 
code word SC as shown in the third line of Figure 3 

The synchronization code word within each batch is followed 
by eight frames of digital data. Each frame is divided into two portions, an 
30 address portion and a message portion. The address code word of the 

message of frame 2 of Figure 3 is shown in line 4 while the message code 
word of the second half of frame 2 is shown in line 5. The address code 
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word is preceded by a digital zero followed by 18 address bus. two (unction 
bits and 10 check bits. The address code word is followed by an even paritv 
bit. The message code word portion of the frame is preceded by a diuital one 
followed by 20 message bits which are followed by 10 check bits and a single 
5 even parity bit. Thus each frame is comprised of 64 bits divided into two 32 
bit sections. 

Synchronization of the base station terminal 200 and the remote 
mobile units 100 is necessary to ensure the units 100 are transmitting at the 
same time that the base stations are listening. Synchronization is also 
necessary to coordinate the division of the large number of remote mobile 
units into groups so that members of one group use different frequency 
hopping patterns from members of other groups. Synchronization of the 
remote mobile units 100 is accomplished by inserting a special frame into the 
POCSAG data which is used to synchronize the units. 

The purpose of synchronization between the reverse paging 
terminal 110 and the remote mobile units 100 is to determine where along the 
pseudo random noise code the frequency hops are to be followed and to 
determine the exact times for transmitting frequencies from within any of the 
hops. This also enables the dynamic changing of a remote mobile units group 
membership such that if one group is experiencing a large number of 
collisions due to simultaneous transmissions, the reverse paging terminal 110 
may re-allocate some of the remote mobile units within that group to new 
groups to minimize collisions. 

Referring once again to Figure 3. eight frames of information 
are transmitted in each burst using the POCSAG format. Remote mobile units 
100 may be assigned to a specific frame within the transmission so that the 
remote mobile units, once recognizing the synchronization code word, can 
scan a specific frame for that remote mobile unifs address. Once the address 
is found, the remote mobile unit can determine any group changes that mav be 
0 required to re-allocate that remote mobile unit to a different group. In 

addition, the POCSAG format is used ,o transmit a fine time synchronization 
code. The fine synchronization code is a transmission of a time pulse at 


20 


25 


an 


WO 96/25673 

PCT/GB96/00270 


12 


exact lime synchronized to a GPS fGlohni iw.- 

* Uilooal Positioning System) clock to 

synchronize all the remote mobile units 100 for r 

iuo lor lime ol transmission. |- or 
example, periodically during the day the r<»v..rc„ 

^ me reverse paging terminal will send a 

synchronization code within the POCSAf: ™a . - 

L5,AO code w ord which is sent at a vcrv 
? precise time. In order to ensure that a nn»ni« . .. 

a precise time pulse is sent, the reverse 
paging terminal 1 10 receives accurate nW ;«r„ 

curate t lmc information using a GPS antenna 

to receive accurate time of dav informitinr. tl 

• iniormat «on. The time used to send the 

synchronization pulse is when the dav ~ i_ 

. ' C, ° ck reaches exactly some multiple of 

0.9 seconds m the preferred embodiment In • - . 

lft , n ... r . f menL In th,s synchronization information. 

IU 20 bits of information are transmitted to oiv^ » 

cu io gi\ e the accurate time of dav 

information. ' 

In each of the remote mobile units 1 00. the microprocessor 
compares this accurate time pulse which will i ndicatc the cxact ^ Qf ^ 

and compare it to its own day clock TK#» „i i 
. . ■ • J The Clock w,thin each microprocessor is 

1 5 accurate down to a few milliseconds, but the time at which the 

synchronization pulse occurs should have a r««i,..- . r- ■ 

nave a resolution much finer than that 

such as down to 0.1 milliseconds for tim*» ^ t .- 

as ror t,me of day. In tins fashion, each of the 

microprocessors in each of .the. remote reverse nanin» a* * 

i reverse paging devices can periodically 

realign its day clock to know within a milli^n^ .u ' 

millisecond the exact time. Each 

20 microprocessor does not actually realinn „u i' l ' . 

y rea,, K n "s clock but changes a clock offset 

withm memory so that it understands how far off its own interna. Cock is and 
can make the adjustment when.using that clock to determine when to start 

transmitting information. 

The synchronization pulse is only transmitted every tew 

25 minutes. However, the resolution of the start n r tu~ ■ 

ui me start oi the message mdicatinii the 

synchronization pulse is very accurate, it being transmitted at 0 090000 
seconds GPS time after a fixed time of day. such as 12:00 GMT. This GPS 
time is accurate to at least within 100 nanoseconds. 

An overview of the transmission format of the remote mobile 
30 unit is shown in Figure 4. The actual transmission of information from the 
remote mobile units 100 is done using Differentia. Bi-Phase Shift Kevin* ■ 
(DBPSK.) modulation on. a frequency hopped carrier of less than one watt 
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The transmission of information from the remote mobile units UK) on the 
frequency hopped carrier may also be done using Frequency Shift Keying 
(1-SK) modulation. The same information is modulated on the two carriers in 
the preferred embodiment of the present invention allowing the base station to 
5 decode the information from either carrier or both carriers. The dual carriers 
are preferably two different narrow hop frequencies. 

Typically a single transmission consists of 53 hops or 53 
changed frequencies selected from a list of narrow band frequencies. The 
frequency selection is based on a pseudo-random noise code list pointing to 
the frequency selection list. The synchronization information tells the remote 
mobile unit 100 where along the pseudo random noise code it should be 
synchronized for transmission of its message and tells exactly the time of day 
so that the remote mobile unit 100 knows exactly when to start transmitting 
the specific frequency so that the base stations BS1. BS2 and BS3. labeled^ 
11 3a, 1 13b and 1 13c. respectively, are looking for that frequency at the same 


10 


15 


time 


20 


25 


50 


In operation. 200 frequencics.are used by the remote mobile . 
unit 100 and the base stations and internally, stored in a list numbered Fl 
through F200. For a specific message. 53 frequencies will be used to transmit 
the entire message. These 53 frequencies are selected based on a 1.000 
member pseudo-random noise code. 

The use of the accurate synchronization signal periodically 
broadcast via the outbound paging signal enables the remote mobile units to 
use lower accuracy components thus reducing the manufacturing cost of 
remote mobile units. For example, high accuracy crystals to track the time of 
day within the microprocessor are available with an accuracy of three parts 
per million. Thus, a time drift of approximately three micro seconds per 
second or 180 microseconds in a minute is the known drift. There are also 
time inaccuracies which are introduced due to variable path length between 
the paging tower to the remote- mobile unit. By employing crystals which are 
cheaper and have an accuracy of the order 50 parts per million, the amount of 
time-of-day drift normally wouldn't be tolerable. However, by using the 
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ic reverse 


synchronization information transmitted on a regular basis from the 
paging terminal, the microprocessor can continually correct its own internal 
day clock so that accurate time of day measurements are always maintained. 
The microprocessor estimates the momentary inaccuracy of the crystal by 
tracking the drift across several synchronization transmissions and dvnamicallv 
adjusts for the frequency drift of the crystal and the offset using internal offset 
registers for accurate time of day information 

Counters are employed within each microprocessor of the 
reverse paging units to compensate for the offset of the frequency based on 
the synchronization time information. There are generally two major factors 
which affect the drift in a crystal: temperature and acceleration. Most of the 
drift is due to temperature, and the remaining drift components are negligible 
The frequency drift in a crystal due to temperature is very slow, on the order 
of 50 Hz over 10 seconds. During a single day the temperature can chanuc bv 
20 or 30 degrees fahrenheit, requiring a time update from the GPS clock/ 
approximately every five minutes. ., 

Remote Mobile Unit Transmit 
The signal sent from the remote. mobile unit 100 to tbe. base 
stations is a spread-spectrum, frequency-hopped transmission using differential 
bi-phase shift keying (DBPSK) modulation on the frequency-hopped carrier to 
transmit digital information. The transmission of information from the remote 
mobile units 100 on the frequency hopped carrier may also be done using 
Frequency Shift Keying (FSK) modulation. The frequency hops are relatively 

slow, the frequencies transmitted are verv narrnu. ,u ♦ 

»^ vcr> narrow and the transmission power 

is extremely small. 

The maximum peak output power of transmission from remote 
mobile unit 100 is limited to less than one Watt to allow use of the 902-928 
MHz ISM band in the United States without the need for licensing the remote 
mobile units as allowed by FCC regulations defined in 47 C.F.R. §15.247. 
Those skilled in the art will readily recognize that other frequency bands and 
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transmissions power levels may be employed depending upon I VC licensinu 
requirements or other frequency licensing requirements of other nationalities. 

The use of an aeeura.e crystal to con.rol each frequency of 
transmission is required within each remote mobile unit 100. For example, 
high accuracy crystals to transmit the narrow bandwidth frequencies used lor 
the frequency hopped transmissions are available with an accuracy of three 
parts per million. At 900 MHz. a 3 ppm drift would place a single frequency 
somewhere within a 2.7 KHz band. To tolerate frequency drift due to aging' 
and temperature, each individual frequency of the frequency hopped signal is 
allocated to a 7.5 KHz band or channel, even though the actual frequency is 
on the order of 200 Hz wide skirt within this 7.5 KHz allocated bandwidth. 
Those skilled in the art will readily recognize that by using alternate 
components, the frequency channels (individual frequency of the frequency 
hopped signals) of 7.5 KHz allocated bandwidth may be wider or more 
narrow depending upon the overall allocated bandwidth for the system. For 
example, 1 KHz or less bands may alternatively be allocated per channel. 

Tests on this system have-shown that by processing the 
received signals at the base stations entirely in the digital domain using the 
combination of unique Fast Fourier Transform algorithms 

to locate and retrieve the frequency hops and by using a 
combination of unique confidence algorithms with a plurality of error 
correction codes, the receiving base station is able to pull the response 
information from a very low power signal from a distance of up to 45 
kilometers (28 miles) in a flat terrain. In a rather noisy urban environment, a 
range of 24 kilometers (15 miles) is the norm. The information within the 
signals is accurately decoded even in. severe multipath and noise conditions. 

As shown in Table 1 . the remote mobile unit message format 
consists, of a preamble and the message body spanning a total of 53 frequency 
hops. Those skilled in the art will readily recognize that longer messages mav 
be transmitted using the preferred embodiment of the present invention, and 
the messages format described here is illustrative but not limiting. Much 
longer message hops to transfer more digital data is also implemented but not 
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described here. Of course, those skilled in the an will readily recognize dial 
shorter messages than those descrihed below are equally possible lor the 
preferred embodiments of the present invention. The message length and 
number of transmission hops arc a matter of desi»n choice. 
5 Thc messa E e preamble consists of a predefined code of ones 

and zeros to gel the attention of the base unit receiver to pull the message out 
of the noise. The preamble consists of 165 bits transmitted across 5 hops, that 
is. transmitted using DBPSK (Differential Bi-Phasc Shift Keying) or 
Frequency Shift Keying (FSK) on five different frequencies selected from the 
10 frequency list with the specific frequencies selected based on the PN (Pseudo- 
random Noise) Code list stored within the remote mobile unit. The sequence 
location within the PN code that the remote mobile unit will begin to follow 
is based on the synchronized time of day. Within a single hop (a single 
carrier frequency), the carrier phase is modulated 33 times to encode the 
15 predefined one-zero pattern of the preamble. 

The message body follows the preamble and consists of three 
groups of data. Each group consists of 30 actual data bits so that the entire 
message is, in the preferred embodiment of the present invention, 90 total data 
bits (although other bit length messages may be chosen). The actual data 
20 encoded within these 90 bits is described above and may be in any convenient 
coded format. Those skilled in the art will readily recognize that a wide 
variety of message formats and encoding of the data bits may be used without 
departing from the scope of the present invention. Thc encoding 

described here, however, has been proven effective, in retrieving the data bits 
25 buried in background noise with a high degree of accuracy and a low actual 
error rate. 

TABLE 1: Remote Mobile Unit Message Fontvat 
Preamble is 165 bits (33 bits x 5 hops) 
Message is 48*33 transmitted bits 
30 (Message is 90 bits actual data) . 

[] = One Frequency Hop 


35 


i — r— i — i — | i i i i — i — i — i — i — i — | — , — . — | — , — , — . — . . . 

preamble I Message Body Spread Over 48 Hops 
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Ouier Mesxnm : Cod i no 
Kach of ihc three groups of message dau. (30 hils each) arc 
BCM1 encoded using a standard 30.63 BCII code and with a single parity bit 
added to form a 64-bit word. This encoding decreases the error rate from 10 - 
5 to I0 ? . This encoding, documented and understood by those skilled in the art. 
can correct up to 6 errors or delect up , 0 13 errors. Detection of corruption of 
a data word thai cannot be reconstructed will cause the base to request a 
second transmission of the acknowledgement message 

10 Inner Coding and lnt,»rUo„j n ,. 

The inner coding of the message will protect the integrity of the 
message with an error rate as high as 25%. Each block of 64 bits of data 
(corresponding to a groups of 30 bits and earlier encoded by a standard 30.63 
BCH code) is split into two sub-blocks of 32 bits (sub-blocks A and B or 
15 Table 2). and a reference bit is added to each sub-block to assist the 
differential encoding to provide a reference bit to the DBPSK or FSK 
decoder. The 33 bit sub-blocks are transmitted over one frequency hop each 
and are replicated 8 times so that the 64-bit block traverses 16 frequency 
hops. In transmission, the 33 bit sub-blocks are interleaved to further reduce 
loss of data, as shown in Table 3, where sub-blocks A and B of Table 2 
correspond to the first group of 30 bits, sub-blocks C and D. correspond to the 
second group of 30 bits. etc. The total message is 53 hops where each hop is 
180 msec in length making the duration of a single message 9.54 seconds. 
TABLE 2 : Interleaving Format for Sub-block 


20 


25 


30 


A = 1 reference bit and 32 data bits = 33 bits 
B = 1 reference bit and 32 data bits = 33 bits 
□ = One Frequency Hop 


B 


ABA 


B 


B 


A B 


ABA 


B A 
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TABLE 3 


10 


15 


20 


Inner Coding and Interleaving of Sub-blocks 

first 


A 
B 
C 
D 
E 
F 


33 bits of is: block 
= second 33 bits of 1st block 
= first 33 bits of 2nd block 
« second 33 bits of 2nd block 
= first 33 bits of 3rd block 
= second 33 bits of 3rd block 

D = One Frequency Hop 


— I — I — I — r 
preamble 
J l l i 


B 


.B 


B 


C 


D 


Those skilled in the art will readily recognize lhat a wide 
variety of data interleaving may be utilized to effect better error tolerance and 
may be substituted for the interleaving described here. Such alternate 
substitute interleaving means are CIRC (Cross Interleaved Reed Solomon 
Code) used in CD (Compact Disc) recording media operating either at the 
block level or at the bit level. 


Simile Hhp F» f ^ nt 
The acknowledgment signals are transmitted by the remote 

25 mobile.. units 1 00 in a 1 .5 MHz band selected from within the 902-928 MHz 
spectrum. The 1.5 Mhz band is divided into 7.5 KHz channels to provide- 200 
channels available in which the frequency hops can occur. Thus, each 
frequency hop is a channel 7.5 KHz wide in which a carrier frequency is 
transmitted. For example, channel one will have a frequency Fl at 

30 902.00375 MHz + -3.75 KHz, channel two will have its center carrier 
frequency at 902.01025 MHz + -3.75 KHz. etc 

Each transmit frequency of each hop will thus be centered at 
the approximate mid-point of the- assigned channel band: however, due to 
inaccuracies in the remote mobile unit circuits and reference crystals, the 

35 actual transmit frequencies will" vary between units. If high quality crystals 

are used to accurately produce the required frequencies, very little drift off the 
center frequency will result. In the preferred embodiment of the present 
invention, low cost crystals are purposely employed to keep the per-unit 
manufacturing costs down. This will allow lor a lower-cost produci sold to 
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the user which will increase market penetration. Thus, reference crystals are 
preferred which have a frequency accuracy of 3 ppm such that at 900 MHz. 
the statistical drill would be approximately 2700 Hz. The crystals center 
frequency within its nominal accuracy also drifts due to aging and temperature 
5 variations, but this drift is slow compared to the transmission times so the 
drift during a single transmission due to these latter variants is unimportant. 

A single frequency hop is shown in Table 4. The 1 5 
millisecond guard time preceding each hop is primarily a settling time for the 
oscillator circuits of the remote mobile units to allow. the internal oscillator 
10 circuit to lock onto the new frequency between hops. Each hop is transmitted 
at a single frequency in which the phase of the carrier is either at 0 degrees 
phase or 180 degrees phase in reference to the phase of the reference bit 
immediately following the quiet or guard time. Thus the first bit is a phase 
reference bit followed by 32 data bits exhibiting either zero phase shift or 180 
degree phase shift to encode the data bits as DBPSK (Differential Bi-Phase 
Shift Keying). In an alternative implementation, each frequency hop may be 
modulated using Frequency Shift Keying (FSK) in which two frequencies are 
used to transmit data bits. One hop frequency may indicate a logical one 
while a second hop frequency may indicate a logical zero. The frequency 
20 shift is minor and the frequency differential is contained within a single hop 
channel. 

Each bit of DBPSK or FSK is a transmission of 5 milliseconds 
of the hop carrier frequency either in phase with the reference bit transmission 
or 1 80 degrees out of phase. 
25 TABLE 4: Single Frequency Hop Format 

Guard time (quiet) = 15 ms 

Single Bit = 5 ms of carrier DBPSK / FSK 

33 Bits plus guard time = 180 ms 

30 


15 


35 


15ms 

5ms 

5ms 

5 ms 

5ms 

Guard 

Ref 

1st 

2nd 

3rd 

Time . 

Bit 

Bit 

Bit. 

Bit 



5ms 

5ms 


31st 

32nd 


Bit 

Bit 
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l-'requcncv Honnim. K.*,y .~.^ 
All of ihc remote mobile units in the market serviced hy the reverse 
paging terminal for message or location llnding use the same pseudo random noise 
code to determine the frequency hops. The p Seu do random noise code is a digital 
5 code of 1000 unique numbers. In the preferred embodiment of the present 

invention, the pseudo random noise code is stored in memory of each of the remote 
mobile units. Those skilled in the art will readily recognize, however, that a linear 
feedback shift register could be used to generate the pseudo random noise code on 
a real-time basis instead of using a look-up table which is presently in the preferred 
10 embodiment. 

The PN (pseudo-random noise) code list is stored in memory and 
maps to a frequency list. In the preferred embodiment of the present invention, the 
PN code list has 1.000 entries which repeat as a sequence. The control means of 
the reverse paging units continuously maintain a count of the proper location within 
15 this list. As described below, the time of day for all remote mobile units in the 
market served by the base terminal are periodically synchronized to ensure 
acknowledgment messages are synchronized to transmit, the hop frequency at the 
proper time and to synchronize the location within the PN code list that each 
remote mobile unit will use to transmit. 
20 The 1.000 member PN code list maps to a 200 member lYcqucncv 

list. In order to allow a large number of remote mobile units to simultaneously 
operate in the same geographic market, the remote mobile units are divided into 
groups and the groups are assigned different sequence segment locations in the 
same 1.000 member PN list. Thus a remote mobile unit from group one will beuin 
25 transmitting a hop at a frequency determined from a ftrst location with the PN 
code, while a remote mobile unit from group two may begin transmitting a hop at. 
a frequency determined from a second location in the PN code. The remote mobile 
units from group one and group two will complete their respective 
acknowledgement messages in 53 hops. Preferably, the sequence of the PN code 
30 used to determine the frequencies of the 53 hops for the remote mobile unit -of the 
first group will not overlap the sequence of the PN code used to determine the 
frequencies of the 53 hops for the remote mobile unit of the second group. More 
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preferable, the frequencies chosen based on ihe non-overlapping segments of the 
PN code list arc orthogonal such that the same frequency is never used bv two 
remote mobile units belonging to different groups. 

In the preferred implementation, the 1.000 member PN code list is 

5 divided into 160 hopping sequences. The remote paging units are divided into 40 
groups with the members of each group synchronized to track the same location in 
the PN code list. The microcontroller of each remote mobile unit, regardless of its 
group membership, continuously runs through the repeating PN code sequence to 
stay in synchronization with the base unit and all other remote mobile units. Each 

10 group of remote mobile units is further divided into four subgroups such that the 
remote mobile units within each subgroup are assigned one sequence within the PN 
code list. Although the 53 hop sequence needed for each acknowledgement 
transmission may overlap the 53 hop sequence used by a remote mobile unit in 
another subgroup, the transmission , sequences of a remote mobile unit of one group 

15 is chosen to not overlap the 53 hop sequence used by a remote mobile unit in 
another group. 

Base Station Desig n 

The analysis of the dual carriers and the decoding of the signals 
20 received by the base stations from the remote mobile, units is done almost entirely 
in the digital domain. The carrier frequencies of the frequency hops are down- 
converted to a lower frequency in each base station and are then digitally sampled. 
The digital samples are then processed to locate the phase information of interest 
for direction finding and message decoding. Each base station is constructed with 

25 a plurality of digital signal processing based receivers which enable simultaneous 
message decoding and direction finding of a plurality of simultaneously 
transmitting remote mobile units. The dual carriers arc analyzed in parallel to 
determine the exact phase separation. Each bases station is also equipped with a 
rubidium atomic clock for accurate time-of-day calculation. Those skilled in the 

30 art will readily recognize that other, less accurate clocks such as GPS clocks can be 
used but the effective location resolution may be affected. 


WO 96/25673 

PCT/GB96/00270 

22 

Differential Time »r Arr; va i Svst«»m 
Figure 6 shows the basic principals of the Differential Time of Arrival 
( DTO A ) System. The example system is configured by at least four receiving 
stations which operate at exactly the same clock. The transmitted signal is received 
5 by the receiving stations at T,. T,. I\. and T 4 respectively. The base stations are 
located at (X,.Y,).(X 2 . Y,).(X,. Y ,) and (X 4 . Y 4 >. The four time 
measurements from the four base stations are resolved in non ambiguous fashion as 
follows: 

10 C a Velocity of Light. 

<1> 

X ,. Y , are the coordinates of the transmitter, and are resolved directly from 
those 3 equations. T, is a dummy variable representing the time of transmission. 

(X.-X.-V + CY ,-Y,)> =c>*(T,-T t )-\ (2) 
15 ( X : - X , ) 2 + ( Y , - Y, ) : = C 2 * ( T , - T ( ) : 

(X 3 -X,) 2 + ( -Y,- Y t ) 2 -C 1 rCT*,-T , t ) = 
(X 4 -X,) 2 + ( Y<-Y t ) 2 =C 2 MT J -T- )^ 

The four equations are dependent (rank 3). Therefore three linear equations could 
20 be from them (by subtracting row I from row 2. row 2 from row 3 and row 3 

from row 4): 

(Define: A P „ = P ,- P, : where P = Oany parameter and CAT ,,= A R, , ) (3) 
2AX 12 X , + 2AY,,'Y, - 2CAR 12 T, - X , 2 + Y,= - X ,= - Y 2 3 - C 2 '(T, 2 -T, 2 )' ' ( 3a) 
25 2 AX 2 , X , + 2AY 53 Y, - 2CAR„ T, = X ,= + Y,- - X ,= . Y,= - C 2 (T, 2 -Tf) (3b) 
2 AX M X , + 2AY„ Y, - 2CAR„ T,. = X :> 2 + Y, 2 - X - Y 4 2 - C 2 (T^ 2 -T 4 2 ) ( 3c ) 

Those three equations could be easily further reduced into two linear equations for 
X , , Y , dependent only upon AR ij. 

30 By dividing equation (3a) by A R,, and equation (3b) by A R : , and 

subtract the results and do the same for equations (3b) and (3c) we get the 
following equations. 
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(Define: D , , ; = X = + Y,- - X , ; - Y;) (4) 
X, J A X |; /A R,, - A X,, /A R,,J + Y, {AY,, /A R,, - AY :; /AR,,', = 

{(D, : - /AR,,) - (D ; r / AR,,) - AR,,)/2 
X, {A X,, /A R,, - A X M /A R, 4 | + V, [AY,, /A R : . - AY, 4 /A R^J = 
5 {(D, ? : /A R, 5 ) - ( D,/ / A R M ) - A R, 4 ]/2 

The location could be resolved with only three bases, but it would 
result in an ambiguous solution which could be resolved by an a priori 
knowledge of the transmitter location. 

10 

Dual Frequency Differenti al Ranp inp Primal.. 
The differential ranging system is based upon transmitting two 
simultaneous frequencies for each hop, with an alternating difference that solves the 
ambiguity. The principal is described in Figure 8. The mathematical derivation of 
15 the differential ranging is based upon the phase difference that is created when two 
frequencies with a small difference in their wavelength are propagating over 
distance between the transmitter and the- receiver. The two signals are represented 
as: ■ ' .. ■ 

20 A, (t) = A„ * Cos (wt + f) (5) 

A, (0 = A,, * Cos ((w-Aw)t + <J> 0 ) (6) 

The phase difference that is propagating through the distance R, 
between the transmitter and the first receiver is: 


A(p, = wR, /C + <j» - (w - Aw) R, /C - <j» „ 


(7) 


The same applies also for the second receiver which is at a distance of R, 
from the transmitter: 

30 

A <p, = wR, IC + $ - ( W . Avv)R; /c . ^ (( (R) 
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Solving Ihc two cqualions: 


AR i: - R, - k, = (A<p, - Acp, )* C/A\v 


(9) 


5 A <p, and A cp ; are measured by the digital receiver at each base station. Aw 

is known a priori. The basic algorithm is as follows: 

1. Compute the phase difference of the two frequencies at the first base 
station. 

2. Compute the phase difference of the two frequencies at the second 
10 base station. 

3. Compute the differential range using formula (9) above. 

Ambiguit y Resolution 
There are two types of ambiguities in the differential ranging system. 

15 First, there is ambiguity that is created if the signal is received by only three 
stations. In that case the hyperbolas are intersecting in two points. 
Mathematically, that ambiguity cannot be resolved unless we have an a priori 
information (for instance - one of the points is in the sea), tracking data (the true 
intersection will move on roads) or other means such as Interferometer Direction 

20 Finding. The second ambiguity is the ambiguity in estimating AR. The reason for 
the ambiguity occurrence in AR is apparent from formula (9) which should be 
expressed as: 

AR12 = R| - Rj = (A(p, - A(p, + 27tn)* C/Aw (jQ) 

25 

where n represents the fact that the difference in wavelength between the two 
frequencies could be larger than 1. Formulating it (10) at the first base station: 

Define M% = (0 - 2*n) 
30 R, = ( n„ + T„ /2tt) * X = ( n l3 + 4',, /2tc) *(?. + A\) t = R, / C 
A,,(t) = A„ * Cos {w( n M + H' u I2n) * X/C.+ <j>) 
\v*>. / C = (w - Aw)*(X + A?.) = 2n 


(12) 
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A,,(i) = A„ + Cos • 2it n,, + M'„ + <t>| 

A,,(t) = A„ * Cos ((w - Aw) t + (()„) = a„ * Cos { 2tt n,, + M' r + <[>„ » 
Acp : = 2* (n„ - n i: ) + M' M - M',, + c|> - «},„ = 2n (n, r n, : ) + AM', + <|> - <}>„ 

5 For the second base station: 

R; = < + H',, /2n) * >. = ( n 3: + M' 3 , /2tc) + AX ): i = R, / C 
A,,(t) = A„ * Cos { 2nn :i + H\, + 4J 

Ay ; (t) = A„ * Cos ((w-Aw)l + <|> n ) = A„ * Cos {2rcn,, + + <)>„ ; 
10 Acp : = 27i(ru,-n, ? ) + H\,-M',, + <t>-<|>o=2n(n,,-n,,) + Af, + 

By subtracting the two measured phase difference we deduce: 

A <p, - A <p : = 2n (n n - n 12 - n 21 + n 2 , ) + AM', - A4'-, = A R,, *Aw/C (1 i) 

15 • ' 

In order for AR ir to be non-ambiguous n„ - - n 2I + n 3: has to be 
0. AR,-,< D. In Figure 7, the length. of one side. of a triangle is longer than the 
; difference between the other two sides. Summarizing it: 

20 Assume Aq>, - A(p, < 7r (14) 

27iAf ♦ AR,v/ C < 7i 
Af < C / (2* AR n ) 
Af < C /(2* D) 

25 For D P raclica »>' limited by 50 Km. Af has to be smaller than 3 

KHz. We chose to use resolution of 2.5 KHz. The problem is that the accuracy of 
AR,, is inversely dependent upon Aw - the larger Aw is. the smaller the error in 
AR,,. 

30 Propagated Wavr Fnm 

The propagated wave form has to resolve two contradicting issues. 
It has a resolution of 2.5 KHz and it has a maximal difference in frequency in 


WO 96/25673 

PCT/GB96/00270 

26 

order to reduce ranging error. The proposed wave form suggests that the ircouencv 
difference Aw will be an integer multiplication of ihe non ambiguous frequency 
S\v: 

5 Aw, = i * ow !lcicN),N = Number or hops. ( , 5) 

The phase difference is incrementing theoretically from hop lo hop by Sep. not 
necessarilv at the same order. 


10 8cp = 8w*AR/C 

Aq>, = AWj * AR / C - i * Sw • AR / C = i*8cp 


(16) 


The evaluation of AR is done by computing ocp which has no 
ambiguities. The proposed detection algorithm is as follows: 

15 1. Create a complex vector ofV,= /V 0 Cos (A(p,).V 0 Sin (Aq»,)J, 

Those values are the quadrate output of the Aq> , calculation and Acp , 
does not have to be calculated explicitly. 

2. Perform FFT ( V k ). . 

3. Find the peak of the FFT and compute its exact "frequency." The 
20 "frequency" is a measure of the phase increments, and provides 

excellent estimation of 87. 

The algorithm automatically improves the signal to noise ratio by the 
factor of the number of hops. There is no dependency between the hop frequency 
and the frequency difference of the number of hops. The process is actually 

25 sampling in the frequency domain. The form of transformation from the frequency 
domain to the time domain is IFFT (Inverse Fast Fourier Transform) rather than 
FFT (Fast Fourier Transform). It has no effect on the results other than 
interpretation. The starting frequency does not have to be 8w or 0 but rather any 
frequency of our choice. It is equivalent to sampling in the time domain which 

30 does not start at T„ = 0 which has no effect on the result of the FFT other than a 
phase shift. 
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Adaptation of Standard NexNtM ™ to nifierential Ranuin«. 
The overall system throughput could be reduced due to two impacts. 
First, the power of each bit is reduced by at least 2. due to the fact that the power 
is divided between 2 frequencies. In Fact the power will be further reduced due to 

5 the non-linearity of the power amplifier which creates spurious frequencies that will 
consume pan of the power. Second, the transmission of two frequencies requires a 
new frequency plan. The frequency resolution is reduced to 2.5 KHz f from 7.5 
KHz) and a special care must be taken to align the hopping with the transmission 
of the two frequencies. 

10 The creation of the two frequencies is done using a regular mixer 

which will multiply the local oscillator by a sine wave in Af/2. The result will be: 

Cos (2«W + * Cos ( 27t{Af/2}t + 4>„ + ) ~ ( , 7) 

0>Cos <2n{f u> + Af/2}t + {^ ln + 4,,,} + ^ 
15 + 0.5*Cqs(27t{f IO - Af/2}l + - Wf . ^ . 

S ±7t/2 (BPSK) 


As is seen, the power is distributed between the two signals. The 
phase modulation is applied through the differential modulating frequency. Both 
20 frequencies are phase modulated in inverse which has no significant since the 
demodulation is differential. The frequency difference is Af. 

Design of the Remote Mohile 1 Init 
In Figure 9. the modulating frequency is generated by the 

25 microprocessor 801 from a port called "TCMP" (time compare) which is capable of 
generating a square wave pulse train. The pulse train is modulated by 0 : or 1 80' 
using a digital XOR gate 802. The modulated square wave is passed through a 
low-pass (or band-pass) filter 803 which clips its harmonics and produces a sine 
wave which produces the two frequencies in the mixer. The ranging transmitter of 

30 Figure 9 is designed to avoid carrier rejection and high order harmonics. The 
power amplifier performance is selected to be extremely linear. This results in 
avoidance of leakage of energy to other channels. 
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The IT circuit for the ranging system is shown in Figure 10. The" 
local oscillator downconvcrts the IF signals from 70 ± 0.65 MHz to 21.4 ± 0.0 1 & 
MIT/, in 36 KHz channels. Two separate converters are used to operate in parallel. 
In order to ensure phase continuity between hops (the phase difference between the 
5 two frequencies must remain between hops) we must base the double loop 
synthesizer on DDS. (Direct Digital Synthesizer). 

The process of finding the phase difference between the two 
frequencies. is based on convolution of the two signals: 

10 Acp = arg( S, <S> S 2 *) (18) 

where arg(*) is defined as the argument of the complex variable (•) The process is 
identical to the process of interferometer DF. The convolution is done on a 
complete hop (45 msec) after the data extraction and the necessary phase reversal. 
15 The phase data of each hop is transferred to the central station. The process of 
differential phase comparison is done in the central station. The location 
processing includes differential ranging computation as well as a combination of 
differential ranging and direction finding. In operation, differential ranging 
requires that periodic calibration be performed in real time: The accuracy of the 
20 system will be in the order of 300 nanoseconds, which mean that liming error in 
the range of 50 nanoseconds must be compensated. The most cost effective way to 
do it is to set up several transmitters in predefined places and to calibrate the 
timing offset accordingly. Assuming a short term drift of 10"'° (high performance 
oscillator), the time between calibrations will be: 
25 T - 50* 10- 9 / 10- ,u = 500 seconds < 19) 

The present system supports remote activation of the calibration 
transmitters. 


30 Accuracy of DTP A 

The accuracy of the DTOA system is determined by comparing it to 
a regular triangulation direction measurement system. The geometric shape 
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representing the points thai have the property of the same difference time of arrival 
relative to a couple of bases is a hyperbola as shown in Figure 7. Accordinu to 
trigonometric law of cosines: 

5 AR = R, -R, = {R- + (D/2)) : + RD Cos (G)» 1 : -{R- +( D /2) 2 - R D Cos(G) ] 1 - 
(2R: + 2(D/2)= - AR 2 ) - = 4 [ ( R : + (D/2)> V - { ( R D Cos (G) >' ] = 

8R- (D/2) : - 4R 2 AR 2 - 4(d/2)- AR : = 8R : (D/2) : - 4 J ( R D Cos (0) j : 
Cos (G) = (AR/D) ( 1 + D ? /4R 2 )'- for (AR/R) « ]. (20) 

10 ln order to simplify the evaluation of the inaccuracy of the location 

due to the inaccuracy of AR estimation, the inaccuracy of the asymptote is 
calculated. This is valid when D/2R«1. Assuming D/2R«1. the asymptote 
formula is: 


15 Cos (G A ) = AR/D 

The asymptote has an ambiguity of ±6 due to its being an even function. 
Estimating the error in G due tothe error in AR is: 


(21) 


20 


o(G A ) = a( AR) / D ♦ Sin(Q A ) (2 2) 


For an example, let E(AR) = 50m. D=20 Km. E(Q A ) is 0.14'. It 
could be shown that the hyperbola error is very close to the asymptote error. The 
error at the transmitter location could be evaluated as: 

25 

a(R) = Ra(G) = o( AR) / {(D/R) Sin(9)} (23) 

Note that the error in the estimation of the DTOA is multiplied by a 
geometric factor (1/ (D/R) Sin(8) ). The constellation of the two receiving bases 
30 and the transmitter could be scaled up or down and the effect on the errors is the 
same. The term (1/ (D/R) Sin (G) ) is recognized as GDOP. 
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C?DOI> e (D/IO Sin(O) (24) 
where GDOP = geometric delusion of properties.- 

Hrror Sour ces in DiftWcntial Ranuin L > 
5 External Thermal Noise is one source of error in the present 

invention. From equation (9) follows equation (25): 

a(AR,, ) = a( A<p,, )* C/Aw (25) 

10 For signal to noise ratio (SNR) above 10 dB the system uses the 

approximation (see Figure 11): 

a(Acp) = 1/ (2 SNR)' - ; ..a ( A<p i: ) = 1/ (SNR) : : (26) 

15 The error in Aq> is dependent upon the algorithm. In order to simplify the 

process, in the alternative, the following is an algorithm for the extraction of the 
frequency out of the FFT. Find the peak of the FFT and : its strongest neighbor.. 
The value of these peaks are marked as P and P'.. The following is the amplitude 
distribution is a sine wave in FFT with boxcar window (signal power equals A„ 

20 which is a function of the signal to noise ratio induced bv Atp ). 

A - A (1 * Sin (71 -8c) . ( 27 } 
N * A<p - 27tk + 5c 

25 

Since Sin (5c) = Sin ( 7t - 5q) the system deduces: 

AVA = 5c / ( 7t - 6q ) (28) 
6c = 7tA7(A + A') 

30 


In order to evaluate the accuracy of the estimation, a point is chosen 
which is the worst regarding the power distribution in the cells of the FFT. At this 
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point A = A\ Assume lhal the noise is while i.e. aA = aA" = \ n . The error in ft; 


is: 


10 


a(5c) = 7i |aA( A 2 + A' 3 )'-/(A + A")-J. = >n/2<2\ * jaA/A] (29) 

A = A „ * Sin ( Tt/2 ) = A „ *2/n 

SNR., = (A ,/ aA) 3 

o( 5c ) = {Ti/2}-/ (2 * SNR„ y- 2 


Since the input noise is distributed over N cells of the FFT- 
SNR W = SNRj * N 


(30) 


SNRj is the input signal to noise ratio. Finally: 
a(AR t2 ) * a( Aq> 12 ) * C/Aw « ( n/8 ) * C / (8f*N w {2 * SNR, }'-) (3i) 
15 where C = 300.000,000 m/s 

N = 53 

5f - 2500 Hz, 

SNRj = 10 dB =10 (integrating over a complete hop) 
20 And the final result of the estimation of AR,-, is: 

a(AR, 2 ) « 27 m ^->) 

Multipath as a Sonrr* r>f P^ r 

25 The effecl of mu, tipath is very different from the effect of the noise 

in two ways. First, the multipath is not necessarily independent, thus the averaging 
that is done through the FFT will not provide an improvement of N 32 . Second, the 
signal to reflection ratio is rarely in the order of 10 dB. Several simulations were 
done in order to evaluate the effect of the rriultipath and the results showed that the 

30 frequency hopping scheme randomizes the effect of the multipath. 
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Dopnlcr S hift as a Source of Frror 
The Dopplcr effect is only on the frequency difference. In 1 Mil/ 


difference we have: 

= v r, = v a iv/" 

(33) 


Af dm .„ = V ina A = V m:ix AI7C 


Where: 

Af drreiN = Maximum Dopplcr shift 

v n. a s = Maximum radial velocity (assume 40 m/s) 

Al* = Maximum frequency difference (up to 1 MHz) 


10 Result in: 


Af dm3X = 0.13 Hz (34) 

Which is negligible since it causes a bias of 0.006 radians in each 
1 5 hop. which results in: 


E{AR nnx }°= E{ A<p J* C/Aw = E{(Ac/N)) * C/Aw = 13.6/nv 
E(p) - Error in parameter estimation 

20 RMU Frequency Inaccuracy as a Int ernal Source of Frrnr 

The modulation of the dual frequency in the Remote Mobile Units 
(RMU) is done with a low cost oscillator which has a possible drift of 100 PPM. 
The error causes a direct error of 0.0001 in the differential ranging (see 32) which 
is an error of several meters in the worst case. Timing is another source of internal 
25 error. The most difficult bias error to overcome is the timing error between 

receivers. It is obvious that even an atomic reference clock with an accuracy of 10 
,: results in 50 nanosecond error within 1 .5 hours. The system design is based on 
calibration transmitters in predetermined places. 

Internal Phase and Frequency Errors is another source of internal 
30 error. Differential frequency errors could be practically eliminated by using a 
TCXO (temperature compensated crystal oscillator) with 1 PPM accuracy. Phase 
errors could not easily be eliminated internally. But. if the frequency generators 
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are phase continuous, the phase error creates -i H.rr, .• . 

. tCS d d,f,er C"l»al ranging bias error (see 

cquat.on I j2 ) which has to calibrate continuously hv the- ,-.n • 
v- , . - - lhc calibration transmitters 

Note the possible group delay at the IF crystal niters. 

* Wave Form 

The transmitted wave form is a H..oi r 
Hiff ( . . . 3 dual fre q u ency transmission with a 

difference o< minimum of 200 KHz to 350 KH7 , - 

^ ^ Hz - w »th resolution of 2.5 KHz. 

The Two-Frequenc y.Differe^^ 

10 The fo "owing description is the com«i«.^ t. 

, K . r complete mathematical analysis of 

the two frequency differential ranging alcorithm t*i • • 

b t algorithm taking ,nto account all known 
sources of error. Let , denote the true time for the /th hop. where , = 0 

corresponds to the beginning of the ;th hop. In the foil™ • 

n u„ . L1 " e follow 'ng equations, a dot 

above the variable indicates its derivative. The tran^min r 

,<-,.. transmitter frequency at the 

15 beginning of the /th hop is given bv 


20 


25 


f*.. - A, +6/ t + 5 f x (i . \ )T 

where . . 

A, is the basic frequency of the /th hop; 
$/ T is the transmitter frequency bias; 
5/ x is the transmitter frequency drift; 
T is the time interval between hops. 

The two instantaneous phases of. the transmitter frequencies are 
♦--^ + /-> + 0.5if.,» + ^ + ^ + |f!W 


(Tl j 


(T2«, 
(T2hj 


where 
30 • 


/.„., is the frequency of the DSB modulation at the /th hop; 

*,« is the initial carrier phase at the /th hop; 

*„,„ is the initial phase of the DSB modulation at the /th hop; 


where 


10 
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», (I) is the phase noise of the transmitter oscillator. 
The instantaneous ranges at the Ah hop are 

«> = + yu> [' + a - \)T\ + o.5« 4M |, + a . 
+ K. + o l0 a -i)7], + 0 .5« t „ 

= ^, + i*.,/ + 0.5a 10 r. 

(T3) 

A' = 1 .2 is the receiver number; 

/•„, is the range between the transmitter and the Ath receiver at the 
beginning of the first hop; 

-V„ is the velochv along LOS between ,„e .ransmiuer and ,he H. receiver 

at the beginning of the first hop: 

«,„ is the acceleration along LOS between the transmitter and the Ath 
I ^ receiver at the beginning of the first hop; 

r,, is the range between the transmitter and the Ath receiver at the 
beginning of the /th hop; 

• v„ is the velocity along LOS between the transmitter and the Ath receiver 

at the beginning of the /'th hop. 

20 The two received signals at the Alh receiver are 

» 0) = 2;« exp [f2n^ h (i . , ^ ( „ + ^ (/f ^ 

25 where 

k„ is the inverse of the speed or light; 
mjij and n hh {l) are the multipaths of the two frequencies; 
ir* u) and ir ttW are additive while noises for the two frequencies. 


0 


In practice, the two frequencies 9r p ~ - j 

quencies are not received separately, but their 

sum ht a) + a) is received. However since the i«,n r™,. • 

«. bince me two lrequencies are separated 

by 7.5 KHz at least, which is much above the bandwidth of each frequency, we 
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assume that ihcy can be separated by prcfillcring. This also explains the factor 2" 
1 \ namely, to relied the fact that the total transmitter power is divided between the 
two frequencies. 

The Alh receiver generates two reference signals whose instantaneous 

5 phases are 

(0 = (fn. + f m .,)U -0 + 6 \ tl ) + (|) rkn . (T5</) 
4Vu (0 = (fiu-L,ja -0 + 6 t rl ) + <J> rU) . (T5b> 

10 where 

0 is the delay, with respect to the beginning of the /th hop. at which 
demodulation starts; 
• 6 t f k is the time bias of the /cth receiver. 

<t> fk 0 is the initial phase of the reference signal of fth -receiver at the rth hop. 
15 Note that the receivers are assumed to have "ideal" frequencies (zero 

frequency deviations), because the reference signals arc synthesized digitally, so/,, 
and /„,, are just numbers in the computer. 

The demodulated signals at the receivers are 

. . 20 " a ^ (t) = ex P (0 ) \. (O (TGa) 

d kh (t) = exp {i2n$ xrh (t)) \ h (i) (T(>h) 

The equations derived so far are used in the simulation. In order to 
derive the estimation let us now ignore the noise. Also express the multipath in 
25 terms of its magnitude and phase, that is. 

>"a„ (0 = A u (() exp {/27i Vjtu (tV (T7a) 

m kb (!) = A kh ft) exp {/2n\|/ w (l)/. (T7h) 

Then, 

30 d ku (0 = T m A ktt (O exp {/2tt [$ rka (t) - ^ (t - Ko r k (t J} - ^n]\ (T8a) 

din 0) = 2 " 12 A * <*> ex P ^ M>n* (i) - 0 - K„r, r/w - * AJ /U];. (T8/» 
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The instantaneous phase of each demodulated siunal contains a 
quadratic term, a linear term, a constant term, and a random term. Denote the 
coefficients of the linear terms by f tlkit and f tllb respectively, li is easy to verify thai 
the linear terms are approximately yiven by 

5 

.U. = k„ \\, (f t l + ./„, , ) - 6/ r 8 \\ (i - 1 )T (T9a) 
Lu, = k« »'o (/:,, - f m , ) - Wr 8 «; 0 - 1 )T. (T9b) 

Suppose we estimate the linear terms (using a Discrete Fourier 
10 Transform with interpolation), and use the estimates to form 

A A A 

<*, =./dl« + f.lih - J.I2a - f.n.h * K«/,..fV,., - V,, ) = K„f„,(Y Xl - V, , ). (T10) 

The use of a , is explained below. 
1 5 Next let us form the signal 

\i (0 = d^Orf'toO) - 0-5A ka (t) A kh (t) 
exp{/27r [<*>,„ 0 - k„ i\ (t)) - <t>, a (I - k„ r t ft)) + <t>, u (l) - (t) + 

* kh (l) -* h .0)\\. (Til) 

20 

We have 


= 2 /„,. *o - %L., 0 - k„ \\, )l + /„,,«,, f t .„ ' : - 2<*> m , (T12) 
25 <fu, W - *w* W = 2/„,, / - 2A„ , & + 2f m , 6/ rf . (T 1 3 ) 

So. 


(l) = Q.5A ttl (OA,* (0 
30 exp {/2ir [2f„,fi,r l .+2f n ,fi„v k j+^^ j . 


(T14) 
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The maximum value of 2f m l K tl v kt t is about 3-10". and the maximum value of ./,„>*„ 
a in r is about 2*10 \ so these two terms can be neglected, leaving 

i k m - 0.5A k ,/t)A u .(Ucxp{j2Tr \2J mt K u r kt - 2 <*W?/^ + ^ (T15) 

» 

5 

Assuming that m ktt (t) and m lh (0 are strongly correlated, wc have i£ Ai , = i£ A/l . Let 
{/„. 0 < « < A' - 1 } be the sampling instants of the ;th hop and form 

(T16) 


p i h = (2n) 'Wangle 


S 5 * ( tr n ) } = 2 f B# s K 0 r ka s -2*^ , -2 f ^ f 6 + 2 f ,6 t rJ , 
ln=0 


10 and 


ft = - A: * -> 2 »-+ 2/^ f (6/ rl - 6/ r2 ). (T17) 

The bandwidth of the signals c/*,/// and d^fij \s limited to about 200 
15 Hz (twice the maximum Dbppler frequency). Therefore they should be low-pass 
nitered to this value and decimated . to a sampling frequency of about 2(KrFtz> 
With the present value of T (50 milliseconds), we get 10 samples per hop. If the 
bit rate is also 10 per second, the range estimation algorithm can be conveniently 
combined with the bit-detection algorithm, since the phase of each bit is needed for 
20 both purposes. The simulation does not include the filtering, but the standard 
deviation of the additive noise is computed in agreement with the assumed 
bandwidth. 

Finally use the estimates a, from equation T10 to compute 


WO 96/25673 


PCT/GB96/00270 


38 


i-l 

/-] 

f=0 

* ~LMt\» - r 2.o) + 2 LA&r x - &t r7 ). (Jig) 

10 As we see. the role of a, is to compensate for the motion during the 

hops and to estimate the initial range difference. Motion compensation depends on 
the receivers carrier frequencies being identical. If they cannot be made identical, 
this compensation will not be performed, and 0, will be used for subsequent 
processing instead of y,. 
1 5 The final stage of the algorithm is the estimation of the range 

difference from {y,. 0 < / < 1 - 1 ) , (or from { p„ 0 < / < 1 - I ) . as discussed 
below). Let t(i) be the ordering of the hops such that f m , (i) = /;,„„ + £ f m „ .0 < f 
< / -1 . Let ti f = exp(/27iy „„) . Then, except for the receiver synchronization 
errors 6 / rk . r| f is given by 

20 

r\ , = cxp{/4« k„ (r, n - r 20 ) (/„„„ + 2 / m „)}. (T19) 

Therefore. 2/„ t0 K (t (r tA) - r 20 ) can be estimated by performing a Discrete Fourier 
Transform with interpolation on the sequence {r| 0 < t < I - 1 }. 
25 If we use (3, instead of y r the sequence of rj f will not be a pure 

sinusoid (even in the ideal case), due to the transmitter motion. However, if f(i) 
, is sufficiently random, the phase deviation seen by (using a Discrete Fourier 
Transform with interpolation) will be random, so this should not have much effect 
on the estimate. For example, if the speed is 30 m/s. the transmitter motion during 
30 the entire sequence of hops is no more than 120 m. With /„„„ + f m „ (/ - 1 ) = 60 
KHz. this translates to about 0.05 •In radians phase error. It is therefore expected 


WO 96A25673 PCI7GB9 6/00270 

39 

that the Discrete Fourier Transform with interpolation procedure will give a good 
estimate of the average range, or the range at the (//2Uh hop. 

The Multipath Model 
5 Multipath is basically a linear phenomenon, that is. the received 

signal is a superposition of several signals, corresponding to several propagation 
paths from the transmitter to the receiver. One of these is the direct path, whose 
impulse response is assumed to be an impulse, with a gain factor a and a zero 
delay (the nominal propagation delay to the receiver is immaterial, and can be 
10 ignored). All the rest are reflected paths. Each reflected path is characterized by a 
random gain and a random delay. These random parameters vary over time in 
general, so the total response at time / to an impulse transmitted at time / 0 is given 
by 

0 

h(t,t 0 ) =ad (t-t 0 ) +£a n (t 0 ) 6 <t-t Q -T r (t 0 )) . {T20) 

15 The case a =? 0 is called Rayleigh fading, While the case a > 0 is called Rice fading: 
In the Rayleigh case there is no direct propagation, only reflections. This happens 
when there is no line of sight between the transmitter and the receiver, and the 
wavelength is too short for diffraction to have an effect. 

When the transmitted signal is s(t) m the received signal is given by 

20 the convolution of s(t) with h(U / (> ). that is. 

r (t) =as(t)+£a r (t)s(t-r n (c) ) . (T2l) 
n 


If we now assume the existence of infinitely many reflected paths, each having an 
infinitesimal contribution. (T21) becomes 


i (t) =as (t) + f°°oc fx, t) s (t-t)dv, 
Jo 


(T22) 


25 where cx(tX) is the infinitesimal contribution at time / of the reflection delayed by 
x. The statistical properties of oefrJ) will be discussed later. 
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Now lcl s(t) he a C\V signal, namely. s(t) =~cxp( \ . Then 


r (c ) - as (C) +f^ a (r * tJexp 1 ititf^i t -r) } ct -cxp ! .7*2* f t .t J [a +m< t ) ; , f rjj) 
where 

m(t) = f°°cc (t, t;exp{-j2iif 0 T} dx . 

As we see from (T24). the effect of multipath on a C\V signal is 
multiplicative, and it depends on the function of a(rj). Here we are interested in 
5 short-term multipath in urban areas and mobile transmitters (or receivers). A 

common model for ct(rj) in such scenarios is as a Gaussian process, stationary and 
narrowband in /. and nonstationary but uncorrected in r. Thus we assume 
EcL(r.t) = 0, 

. r(x, , x 3 ; /, , / 2 ) = Eoc(r } J x )a(x 2 , /,) = r, (*, - / 2 ) T t (x x )bix t - x 2 ). (T25) 
10 The temporal correlation T, (7, - / 2 ) depends mainly on the physical 

motion of the transmitter. In particular, if the doppler frequency of the transmitter 
is f tlttpplcr (corresponding to the wavelength and the transmitter speed), then the 
spectral density of T, (7, ~ h} is limited to ±f tUtpptcr . The shape of the spectral 
density is difficult to estimate, since it depends on the given scenario: mainly on 
15 the spatial distribution of the reflections. In the simulation we used a flat spectrum 
in the range ±f dopphr for simplicity and "universality." 

The function T t (x,) describes the relative multipath power at 
different delays. It is also scenario-dependent and difficult to estimate. However, 
based on the graphs known in the art, we assume a functional dependence of the 
20 form 

Mx^exp'-T./T^!, (T26) 
where t is a scenario-dependent parameter. Note that when T T fr,) is plotted in 
dB. it shows a decreasing straight line. 

It is interesting to compute the multipath correlation for both 
25 different times and different frequencies. Let m(/, , /,) be the multipath at 
frequency /, and time /, . and similarly for m(t 2 . f 2 ). Then 
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r. A t w t,; f, t f-,) =Em{ t, , f.-)/n* ( t :: , 


E S°° IT" tj) a (l2 ' t;?) eXp * "J 211 < f i t l"^ X 2) 1 cfXjdx- 


= rf° I\ ( t 2 -t 2 )V x (x,) 6 (x 1 -t 2 )exp{-j2n ( f,x, -f 2 x 2 ) } dx.dt- 
J o J c 


=j~r t ( .t J -cpr t <t J )exp{-j2ii (f,-/,) left, 
10 0 


Fc(tj ^IsE—... (T27) 


1+J271 (f J -f 2 >T 


JTip 
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Table 1: Simulation Results lor the Two l ; rcuucncv Model 


/; 

SNK 

Ml' 

SI 1) 

Min 

Max 

60 MHz 

6 dB 

0 

8 m 

130 m 

633 in 

60 MHz 

9 dB 

0 

-2 m 

123 m 

561 m 

60 MHz 

6 dB 

1 

2 m 

80 m 

328 m 

20 MHz 

6 dB 

0 

-56 m 

334 m 

1530 m 

20 MHz 

9 dB 

0 

4 m 

211 m 

801 m 

20 MHz 

12 dB 

0 

4 m 

152 m 

344 m 

20 MHz 

6 dB 

1 

6 m 

208 m 

921 m 

20 MHz 

9 dB 

1 

-21 m 

140 m 

831 m 


Table 1 shows the results of simulation of the present system with 
the given parameters. The frequency is selected as described above as/„,„ + 
f m0 (/ - 1).. SNR is the signal to noise ratio. MP = 0 indicates Rayleigh fading 
15 and MP = 1 indicates Rice Fading. STD is the standard deviation (one a), and the 
resolving power of the algorithm is given us minimums (Min) and maximums 
(Max). 

20 Although specific embodiments have been illustrated and described herein, 

it will be appreciated by those of ordinary skill in the art that any arrangement 
which is calculated to achieve the same purpose may be substituted for the specific 
embodiment shown. This application is intended to cover any adaptations or 
variations of the present invention. 
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CLAIMS 

I. A mobile transmitter location system, cotnprisint* : 

a mobile transmitter having a unique identifier, ami having a 

transmitter lor simultaneously sending a first frequency-hopped spread spectrum 
5 carner^ignal and a second frequency-hopped spread spectrum-carrier siunal, at 

least one of which containing the unique identifier, the first frequency-hopped 

spread spectrum carrier signal having a different frequency than the second 

frequency-hopped spread spectrum carrier signal; 

a first base station having an antenna operable for receiving the first 
10 and second frequency-hopped spread spectrum carrier signals and capable of 

determining a first time of arrival of the frequency-hopped spread spectrum carrier 

signals; 

a second base station also having an antenna operable for receiving 
the first and second frequency-hopped spread spectrum carrier signals and capable 
15 of determining a second time of arrival of the frequency-hopped spread spectrum 
carrier signals; and 

a central station connected to the first base station and the second 
base station, and operable for: 

.1.) receiving the first time of arrival of the first and second 
20 frequency-hopped spread spectrum carrier signals from the first base 

station; 

2.) receiving the second time of arrival of the first and 
second frequency-hopped spread spectrum carrier signals from the 
second base station; 

25 3.) generating a first hyperbolic plot based on the differential 

time of arrival of the first and second frequency-hopped spread 
spectrum carrier signals at the first and second base stations; and 

4.) using a third information source to determine the location 
of the transmitter along the first hyperbolic plot. 

30 

2. The mobile location system according to claim 1 wherein a reference station 
operates periodically to transmit a reference first and second frequency-hopped 


44 

spread spectrum carrier signal from a known location 10 calibrate the first and 
second base stations and the central station. 


3. The mobile location system according to claim 1 wherein the first base 

5 station and the second base station both include atomic clocks to synchronize the 
time of arrival calculations. 

4. The mobile location system according to claim 1 wherein the third 
information source is a geographic map of the area around the base stations which 

1 0 provides information on the probable location of the transmitter along the first 
hyperbolic plot. 

5. The mobile location system according to claim 1 wherein the third 
information source is a third base station which calculates the time of arrival of the 

15 first and second frequency-hopped spread spectrum carrier signals and wherein the 
base station generates a second hyperbolic plot to intersect the first hyperbolic plot 
to identify the. location of the transmitter. 

6. The mobile location system according to claim 1 wherein the transmitter 
20 transmits a plurality of hops each having a different first and second frequency- 
hopped spread spectrum carrier signals and the time of arrival is determined by 

a.) calculating the phase difference between the first and second frequency- 
hopped spread spectrum carrier signals for each of the plurality of hops to produce 
a plurality of phase differentials: and 
25 b.) performing an Inverse Fourier Transform on the plurality of phase 

differentials to produces a time of arrival. 

7 - A method of locating a mobile transmitter location, comprising the 

steps of: 

30 receiving at a first base station a dual radio carrier having two 

frequencies and calculating therefrom a first time of arrival of the dual carrier: 
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receiving at a second base station a dual radio 
carrier having two frequencies and calculating therefrom 
a second time of arrival of the dual carrier; 

comparing the first time of arrival with the second 
5 time of arrival; and 

determining the location of the origin of the dual 
radio carrier by hyperbolic intersection location. 

8. A mobile transmitter having a unique identifier for 
10 use in the method of claim 7 comprising a transmitter for 
simultaneously sending a first frequency-hopped spread 
spectrum carrier signal and a second frequency-hopped 
spread spectrum carrier signal, at least one of which 
containing the unique identifier, the first frequency- 
15 hopped spread spectrum carrier signal haying a different 
frequency than the second frequency-hopped spread 
spectrum carrier signal. 

9 . A base station for use in the method of claim 7 
20 comprising an antenna for receiving the dual radio 

carrier and means for determining the first or second 
time of arrival of the dual carrier. 

10. A mobile transmitter location system comprising: 
25 a first base station having an antenna for receiving 

a dual radio carrier having two frequencies and means to 
calculate therefrom a first time of arrival of the dual 
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carrier; 

a second base station having an antenna for 
receiving a dual radio carrier having two frequencies and 
means to calculate therefrom a second time of arrival of 
5 the dual carrier; 

comparing means for comparing the first time of 
arrival with the second time of arrival; and 

location determination means for determining the 
location of the dual radio carrier by hyperbolic 
10 intersection location. 
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